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>y 1 . Introduction The P-L relationship for LPVs has been investigated us- 

^, . , T s , • ing the Hipparcos output catalogue ( Ivan Leeuwen et alJI 19971) 

This paper investigates the period-luminosity (P-L) relation- ^ , ^ u- . . j- » ti, t 

^ ^ , f ^ . . ror samples or obie cts at the same d istance (in the Large 
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Abstract. We analyze the K band luminosities of a sample of galactic long-period variables using parallaxes measured 
by the Hipparcos mission. The parallaxes are in most cases re-computed from the Hipparcos Intermediate Astrometric Data 
using improved astrometric fits and chromaticity corrections. The K band magnitudes are taken from the literature and from 
measurements by COBE, and are corrected for interstellar and circumstellar extinction. 

The sample contains stars of several spectral types: M, S and C, and of several variability classes: Mira, semiregular SRa, and 
SRb. We find that the distribution of stars in the period-luminosity plane is independent of circumstellar chemistry, but that 
the different variability types have different P-L distributions. Both the Mira variables and the SRb variables have reasonably 
well-defined period-luminosity relationships, but with very different slopes. The SRa variables are distributed between the two 
classes, suggesting that they are a mixture of Miras and SRb, rather than a separate class of stars. New period-luminosity 
relationships are derived based on our revised Hip parcos parallaxes. The Miras show a similar period-luminosity relationship 
to that found for Large Magellanic Cloud Miras bv lFeast et alj il989l) . 

The maximum absolute K magnitude of the sample is about —8.2 for both Miras and semi-regular stars, only a little fainter 
than the expected AGB limit. We show that the stars with the longest periods (P > 400 d) have high mass loss rates and are 
almost all Mira variables. 

Key words. stars:carbon - stars:AGB and post-AGB - stars :pulsation 



ship, or relationships, for long-period variable stars (LPV) in Mage llanic Cloud jOlass & Llovd Evans. ..198lUFeast et alJ 



the solai- neighborhood. P-L relationships for stars in the in- iggQ^ lGroenewegen & Whitelockir996HCioni et alJ200ll) .7;;d 
stability sti-ips of the HR diagram have proven to be of fun- the bulge of the Galaxy jAJMet al.-2001)). The last two pa- 

damental value in using variable stars as secondary distance , j, ^, ^ j ^ w 1.1. 

■' pers make use of the vast data compilation on variable stars 
indicators. Stellar variability is also a powerful tool for study- ■^ c ■^ ■ t-i. ^ j n t 

■> ^ J available from microlensing surveys. These papers nnd P-L re- 

ing stellar structure, chemistry and mass loss (e.g. Wood 1990 : , ^. . • ^ Tl^/t/-^^/r• -^uu *u j u u 

I — — II -\ I • 11 ^ I ^ ° — ^ lationships for LMC Miras with both oxygen and carbon chem- 

IWillsonll200(]t lAlvarez et alJI2001l) . The present paper investi- • ^ j w u- f 1 • ui u- u 

' ^' ^ ^ ^ ^ istry, and relationships for semi-regular (SR) variables which 

gates whether one or more P-L relationships exist for LPVs in , , f iv^- u ^ ti, » 

° ^ have the same slope as for Miras but with two sequences sug- 

the local region of the Galaxy, using re-derived Hipparcos par- ^ ,^ j^.ui ju. j uur 

o ff f gested to correspond to t he long a nd short periods which of- 

allaxes, A magnitude data from the literature, and A -band flux ^ • . • . J„. , ,^nr\^^ ^ j- 

, ^ , ^ ten coexist in these stars ICioni et aLc2001.): see also th e dis- 



densities and variability data derived from the COBE DIRBE Hussion bvlSzatmarv et all JT994 . bedding et alJ 1199^), and 
time-series data ( Boggess et al,1992, . Hauser et al,1998,) . Izijlstta et al.V2002')] . In the Galactic Bulge jAlardetaylool 
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their effective temperatures low, < 2500 K, so that the bulk of 
their flux is emitted at near-infrared wavelengths, where inter- 
stellar extinction is much lower than at optical wavelengths) 
the establishment of a well-defined P-L relationship is of sig- 
nificant interest. Further, the P-L relationship establishes the 
mode in which the star is pulsating, and gives information on 
the effect of pulsation on mass loss. 

Measurements of the P-L relationship in the LMC and the 
Galactic Bulge provide the slope a of the relationship: 

= alogP(days) + 6 (1) 

and the zero-point b (using the kno wn distances to the LMC 
and the bulge). However, models (e.g. IWoodI 1 99(1 12000 ) show 
that the zero-point absolute luminosity is metallicity depen- 
dent. Comparisons of the P-L relationship derived for nearby 
galactic LPVs, whose absolute distances are known from par- 
allax measurements and for which there is much more infor- 
mation available on metallicity and mass loss, are thus worth 
investigating. 

Establishment of P-L relationships for nearby LPVs is dif- 
ficult, however. Even after Hipparcos ( ESA 199 7) there are 
relatively few well determined parallaxes for these rare, lumi- 
nous starSj_and their use requires attention to statistical biases 
(cf. lArenou & L uri 1995). Several studies have assumed that 
the P-L relationship for ne arby Miras has the same slope as 
found for the LMC Miras (Groeneweaen & Whitelock 1996; 
Bedd ing & Ziilstra 199 8: Whitelock et al 2000; Bergeatet al,, 
200 ll) and have used this relationship to calcul ate the zero point 
for galactic populations. 'Barthes et alJ ( Il9 99^. however, using 
both direct and statistical parallaxes for a sample of Miras. find 
a diffe rent slope than that seen in the LMC, while lAlvarez et af] 
11997^ identify two groups of Miras with different zero points 
(cf. also Jura & Kleinmann 1992). 

Hipparcos provided the first parallaxes for most LPVs for 
which this information is now available. However, most of the 
parallaxes available even from Hipparcos are measured at low 
confidence. Added to this are reduction problems associated 
with the very red colors and extreme variability of these ob- 
jects. These factors make the evaluation of absolute magnitudes 
for these objects quite uncertain. 

We have recently undertaken a re-analysis of the Hipparcos 
astrometric data for red giant stars with the aim of improv- 
ing the accuracy of the parallax information for t hese objects. 
The re -reduction scheme is described in detail bv'Platai s et alJ 
1I2OO3I) and Pourbaix et al.. (2003 ). The present paper uses im- 
proved parallax information from these studies, plus some 
Hipparcos catalogue parallaxes, to re-examine the absolute 
magnitudes of variable AGB stars. We work entirely in K 
magnitudes, using magnitude data from the literature and from 
COBE-DIRBE data. We also incorporate corrections for inter- 
stellar and circumstellar extinction, and use the DIRBE data 
to examine the K band variability of the stars. The assem- 
bly of the data set is described in the next Section. Section 
3 discusses the distribution of LPVs in the P-L plane and in- 
vestigates their P-L relationships. We also examine the maxi- 
mum AGB luminosity for both carbon- and oxygen-rich stars. 
Section 4 briefly discusses data on period variabiHty, from the 



literature and from the DIRBE data, and discusses the relation- 
ships among the different types of LPVs. Section 5 discusses 
Mira variables with long periods and large mass loss rates. 

2. Data 

Table 1. Stars for which a new solution based on the Hipparcos 
observations was derived, w and £{vj) are the parallax and 
its error in mas. w — whip is the difference between the re- 
vised parallax and the published one, Ak is the total interstellar 
and circumstellar extinction in the K band (see text). K is the 
2.2/im apparent K magnitude. P is the period in days, 'var' is 
the variability type from CGVS, 'Chem' is the carbon/oxygen 
chemistry of the atmosphere 

HIP GCVS w e{w) -uo - ronip Ak K P var. Chem" 

The sample of objects discussed herein was basically de- 
fined by the available data: we need parallaxes, K magni- 
tudes, periods and variability types. Secondary information in- 
cludes the position (which together with the parallax allows 
the calcul ation of interstel lar extinction); the 12/im flux den- 
sity ( Beichman et aljfl988l) which, together with the 2/im flux 
density, gives an estimate for the circumstellar extinction, and 
the circumstellar chemistry, since the dust opacity at 12/im is 
strongly dependent on its composition. 

The data are given in Tables [2 and |2] where we list: the 
Hipparcos catalogue number; the variable star name; the paral- 
lax w in milliarcseconds (mas); the parallax error £(117) in mas; 
the total (interstellar and circumstellar) K band extinction: the 
K magnitude; the period in days; the variability type; and the 
chemistry (O = "oxygen-rich", i.e. n(0)> ?i(C); C = carbon 
star and S = S star). The sources of these various quantities are 
discussed below. 

2. 1. Parallaxes 

The sample of LPVs analyzed in this wor k was selected from 
stars observed by Hipparcos ilRSAlll997h . This mission pro- 
vided accurate astrometric measurements for a sample of some 
120000 stars which is approximately complete to a brightness 
limit of V ^ 8 but also contains many fainter stars. Included 
in the Hipparcos input catalogue were several hundred evolved 
cool giant stars, and Hipparcos provided the first measured par- 
allaxes for almost all of these. However, because of their rela- 
tively large distances, variability, red colors and, for some stars, 
their faint apparent magnitudes, the parallax errors for these 
stars are often larger than the nominal Hipparcos 1 a accuracy 
of 1 mas. We have therefore undertaken a re-analysis of the 
Hipparcos Intermediate Astrometric Data (IAD) for these ob- 
jects. The re-processing scheme contains two components: de- 
riving epoch V — I p hotometry for each IAD using the method 
described by IPlatais et al . (2003). and then reprocessing the 
IAD to derive parallaxes using epoch chromaticity corrections. 
This accounts for the color dependence of the Hipparcos detec- 
tor, a particularly important aspect of the data processing for 
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Table 2. Same as Table [T]but for SRb variables for which the 
Hipparcos data have not been reprocessed. 
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RR Eri 
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1.04 
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0.52 


97 
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CY Eri 
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0.01 


1.96 


25 
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RV Cam 


3.01 
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0.06 


0.52 
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W Ori 
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0.01 
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0.02 


-0.46 
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34922 


LZ Pup 


16.46 


1.27 


0.02 


-2.49 


140.6 
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RC CMi 


7.77 


0.99 


0.00 


0.83 


35 
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RT Hya 


3.68 
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0.01 


0.09 


290 





42502 


AK Hva 


5.12 


0.92 


0.02 


-0.61 


75 





44050 


RT Cnc 


2.94 


1.11 


0.01 


0.20 


60 
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55639 


T Crt 


5.04 


1.06 


0.01 


1.20 


70 
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BK Vir 


5.68 


1.12 


0.01 


-0.88 
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62611 


SV Crv 


2.72 


0.97 


0.01 


1.07 


70 
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RT Vir 


7.25 
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0.01 


-1.12 
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64569 


SW Vir 
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0.01 


-1.87 
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0.76 


0.02 


-0.72 


90 





71644 


RV Boo 


2.54 


0.98 


0.00 


-0.01 
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2.97 


1.19 
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0.12 
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7.26 
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0.01 


0.67 
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XX Dra 


3.52 


0.56 


0.01 


1.38 
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82249 


AH Dra 


3.39 


0.7 


0.01 


0.57 
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87850 


OP Her 


3.26 


0.54 


0.02 


0.05 


120.5 





89669 


IQ Her 


3.82 


0.87 


0.04 


0.04 


75 
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91781 


V3879 Sgr 


2.67 


1.01 


0.05 


-0.42 


50 





92862 


RLyr 


9.33 


0.52 


0.01 


-2.10 


46 





96204 


V450 Aql 


3.8 


0.9 


0.04 


-0.08 


64.2 





97151 


V973 Cyg 


3.6 


0.57 


0.02 


1.49 


40 





100935 
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3.44 


1.19 


0.01 


-1.59 


347 
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5.28 


0.63 


0.02 


-1.48 


131.1 





107487 AG Cap 


2.64 


0.85 


0.01 


1.20 


25 





110099 UWPeg 


3.38 


1.57 


0.01 


2.27 


106 





118249 


SPhe 


3.93 


0.87 


0.03 


1.56 


141 






the very red, and highly variable, LPVs. For the present study, 
we used three data sets based on these new reductions: carbon 
stars - the requirement of a measured period remo ves the much 
less luminous R type stars from the sample (cf. iKnapp et alJ 
I2OOI1) . leaving the cooler, more luminous, variable N-type car- 
bon s tars; Variability-Induced Movers (VIMs; Pourbaix et al. 
I2OO3I) : and oxygen-rich Miras. The parallaxes and their errors 
are given in Tabled Note that in many cases the parallax is not 
significant. 

The samples of carbon stars, VIMs and oxygen-rich Miras 
defined above are complete in the sense that they comprise all 
such objects observed by Hipparcos except for the small num- 
ber without measured periods or K magnitudes. The Hipparcos 
carbon star and VIM samples contribute a fair number of SR 
variables: to round out the data sets, we extracted from the 
Hipparcos catalogue data for the remaining oxygen-rich SRb 
observed by Hipparcos. The data for this sample, with the same 
information as in Tabled are given in Table 2. 



0.15 - 




1000 



Period (days) 

Fig. 1. Frequency of periods from the CGCVS. Heavy solid 
line: Mira. Light solid line: SRa. Heavy dotted line: SRb. Light 
dotted hne: SR. 



2.2. Periods and Variability Types 

The monumental Combined General Catalogue of Variable 
Stars (CGCVS, Kholopov et al. 1998) contains data for many 
types of variable stars, including amplitudes at visual wave- 
lengths, periods, and classification. LPVs are classified into 
several subtypes: Miras, with long, well-established, repeat- 
able periods (typically 300 d or longer) and large amplitudes 
(AV > 2.5 magnitudes): semiregular variables with shorter, 
less repeatable periods and smaller amplitudes, which are sub- 
classified into types SRa and SRb: and irregular variables, type 
Lb, which show no well-defined period (in fact, periods have 
been measured for a few of these, and thei r variability behav- 
ior is indistinguishable from that of SRb; IKerschbaum et alJ 
I2OOII) . Many semiregular variables show two or more periods, 
and some switch between two or more periods. We discuss 
these variability types further in Section 5. 

Figure 1 shows the period distribution for LPVs from the 
CGCVS. Many of the semiregular variables are classified only 
as SR. The period distributions for all types cover a similar 
range, and there are few stars of any type with periods longer 
than 400 days (we discuss this further in Section 5). The me- 
dian period of the Miras is more than twice that of the SRb. 
The stars classified only as SR have a period distribution clos- 
est to the SRb. The SRa distribution is intermediate between 
the Miras and SRb. In this paper, we will use three classes: 
Mira, SRa and SRb, combining the last with the SR and the 
very small number of stars classified as Lb for which a period 
has been found. 

The periods and variability types in Table [0 were taken 
from the CGCVS and the Hipparcos variabihty annex. Periods 
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Fig. 2. K magnitude from DIRBE 2.2 /im flux densities vs. K 
magnitude from the literature. The dashed line shows equality. 
The bars on the DIRBE magnitudes show the observed range 
of magnitude variation over the time span of the COBE obser- 
vations. 



measured by Hipparcos agree well with the CGCVS periods in 
almost all cases. 



2.3. K magnitudes 

We analyzed the P-L relation for LPVs using K magnitudes, 
which are available for almost all bright LPVs from the lit- 
erature. The variation of Miras between maximum and min- 
imum light can exceed 10 magnitudes at visual wavelengths, 
since the visual luminosity is a highly non-linear function of 
photospheric temperature because of molecule formation (for a 
recent discussion see ,Reid&_G oldston 2002); while the vari- 
abiht y amplitude is much smal ler (about 0.5 magnitudes) at K 
(e.g. IWhitelock & Feast "2000V The K luminosity is also less 
affected by chemical composition ( Wood 2(X)0) and is near the 
peak of the stellar spectral energy distribution. 

K magnitudes were taken from the literature, includ- 
ing the IRC (,Neugebauer&Leighton 1969.). the com pila- 
tion bv lOezari et all (Il999l) . iKerschbaum & HrorJ ( 1 1994 and 
IWhitelock & Feast' ('2000V Except for the last reference, which 
gives K magnitudes averaged over the light curve, these are 
single-epoch observations. 

To study the variance of the K magnitudes caused by 
variability, we extracted 2.2/im light curves from the data 
from the Diffuse Infrared Background Explorer (DIRBE - 
[Hauser et al.'l998') on the Cosmic Background Explorer COBE 
(Bog gesseLal . 1992). DIRBE provides well-sampled partial 
light curves (sometimes with several observations per day) for 
bright point sources at wavelengths from 1.25/im to 240/im 
over a time interval of almost 300 days (11 December 1989 to 



100 200 
At (days) 



300 



2.5 



E 
3, 



1 .5 - 



1 - 



rPo° o 



100 



200 

Period ( 



300 400 
days) 



500 



Fig. 3. Ratio of maximum and minimum 2.2 /im flux densities 
versus (a) time between maximum and minimum 2.2 light, 
as measured by DIRBE (upper panel) and (b) period from the 
CGCVS (lower panel). Filled circles: Mira variables. Open cir- 
cles: SRb and SR variables. Crosses: SRa variables. 



21 September 1990). The gaps in the data and the short du- 
ration of the COBE mission relative to the typical period of 
an LPV limit the usefulness of COBE for measuring periods 
(see e.g. Smith et al. 2002). Because of the large DIRBE beam 
(40') the data are useful only for the brightest point sources, 
those brighter than about 0-1 magnitudes. 

The DIRBE 2.2/im light curves were extracted for 
the brightest objects in Table ^ (approximately mag 
and brighter at K) in the present sample using software 
written by N. Odegard and D. Leisawitz and available at 
http://cobe.gsfc.nasa.gov/cio/browser.html 
The light curves of some of these stars will be discussed in 
Section 4. Figure 2 shows the DIRBE K magnitudes versus 
those from the Uterature. The DIRBE 2.2/im magnitudes 
were calculated from the observed flux densities assuming a 
value for the flux density of a r nag star of 654 Jy [usi ng the 
calibrations of a 0.03 mag star bv lBessell & Bretll ( Il988h l. 

The DIRBE time series data were used to calculate an aver- 
age flux density for the entire data set plus the mean deviation 
from that average. In most cases most of the deviation is found 
to be due to stellar variability rather than noise. The ranges in 
magnitudes for the COBE data are also shown in Figure 2. 
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In general the agreement between the data from the liter- 
ature and from DIRBE is excellent. The mean internal devia- 
tion from the COBE data is 0.2-0.4 magnitudes for Mira vari- 
ables and about 0.02 - 0.1 mag for SR and SRb variables. In 
almost all cases the difference between the DIRBE and the lit- 
erature magnitude is within these deviations, which therefore 
dominate the measurement noise. The data for a few stars show 
larger discrepancies, and may be due to longer-term variations 
in these stars. 

Figure 3 shows the ratio of the maximum to minimum 2/^im 
flux density for stars for which a sufficient portion of the light 
curve is observed by DIRBE to observe at least one maximum 
and one minimum (only a minority of the stars in Table 0. 
Figure 3a shows this ratio versus the time difference between 
the minimum and maximum, or vice-versa - not necessarily 
half the period, since the light curves for LPVs are often asym- 
metric. Figure 3b shows the ratio versus the period from the 
CGCVS. Figure 3 shows that Miras are far more readily dis- 
tinguished from SRb by amplitude than by period (cf. also 
ICioni et al. 2001). although stars with (tmax ~ ^min) shorter 
than about 70 days are always SRb. The SRa variables overlap 
both the Miras and SRb and are therefore likely to be a mixture 
of the two types (cf. Figure 1). 

In the next section, then, we discuss only two samples: 
Miras and SR/SRb. We do not include the SRa in the analysis. 
We assume that the error in the COBE median K magnitude 
is 0.1 mag for SR and 0.3 mag for Miras,in both cases dom- 
inated by variability; these are the median values found from 
the COBE light curves. 

2.4. Extinction 

The extinction towards the stars has two com ponents, inter- 
stellar and circumstellar. We used the model of iHakkila et alJ 
ill 9971) to calculate Ay, the total interstellar extinction in V, 
from the position of each star and its distance, estimated from 
w~^. The K ba nd extinction is then found from Ak^ = 
0.114 Av (Carde lli et al.ll989l) . 

Many of the stars in TableHja re losing mass and surro unded 
by dusty circumstellar envelopes. Ivezic & Elitzur ( 1997) show 
that the resulting spectral energy distribution is a function of 
the dust composition and the optical depth. The circumstellar 
extinctions for the stars discussed in this paper are not large 
by the standards of mass-losing AGB stars (they are all visible 
objects) and the optical depth can be estimated from the ratio 
of the 12/im flux density (emitted almost entirely by the dusty 
envelope) to the 2fim flux density (mostly emitted by the pho- 
tosphere). The relationship between the 12/im/2/^m flux ratio 
and the extinction at K was estimated using a spherically sym- 
metric dusty envelope model ( Knapp et al 1993) using sil icate 

grains for oxygen and S stars fcf. llorissen & Knappll998l) and 
carbonaceous grains for carbon stars. The model assumed con- 
densation temperatures of 800 K an d 1000 K for silicate a nd 
carbonaceous grains respectively (cf. Ilvezic & ElitzuJ 19971) to 
define the radius of the inner edge of the model envelope. The 
model's central star is a blackbody with — 2500 K, = 



5 

5000 L0, for which the intrinsic ratio is S12/S2.2 = 0.18. For 
low optical depths 

Aif, « 0.07 (^-^-0.18^ (2) 
for silicate envelopes, and 

Ak^ « 1.1 - 0.18^ (3) 

for carbonaceous dust. The large difference between the two is 
because of the strong contribution of the silicate 9.7/^m emis- 
sion feature in the IRAS 12/xm band. A very small number of 
stars have no IRAS observations: we assume a circumstellar 
extinction of mag for these stars. 

The total (interstellar -1- circumstellar) extinction at K band 
is listed in Tables ^ and |2l The uncertainties in both the inter- 
stellar and circumstellar extinction are large, but the contribu- 
tion of these errors to the errors in the computed absolute K 
magnitude are dwarfed by those in the parallax and by stellar 
variability, and thus will be neglected in the remainder of this 
paper. 

3. The Period-Luminosity relationship 

The stars listed in Tables Handle often have poorly measured 
parallaxes. Figure |4] shows the distribution in the period- A/^f 
plane of stars with reliable m, defined by: 

- probability that n7 > at the 95% confidence level. The 
mean Hipparcos parallax error is ^ 1 mas with some small 
scatter around this value depending on the location of the 
star on the sky and on its brightness. LPVs are very lumi- 
nous, so many observed by Hipparcos may be at distances 
> 1 kpc. The probability that the parallax is detected is de- 
termined by including/not including the parallax in the as- 
trometric solution, and defining the probability that the par- 
allax is detected using the F-test (cf. Pourbaix et al.,2003h . 

- n7/e{Tu) > 2 

- nj > 0. 

- e{-nj) < 2 mas. As discussed bv lPlatais et al.l ll2003l) the er- 
rors on the Hipparcos parallaxes for some relatively bright 
red giant stars are anomalously large, and measurements 
with large errors often produce inaccurate parallaxes. 

The distance is calculated from D — lOOO/nj and the ab- 
solute magnitude Mk from this distance and the observed K 
magnitude, corrected for circumstellar and interstellar extinc- 
tion. The errors on Mk are assumed to be due to the errors 
(or rather range of photometric variation) in K and the paral- 
lax error. In |4] the stars are coded by variability type (Mira, 
SR/SRb and SRa) and by circumstellar chemistry (combin- 
ing the oxygen and S stars). Also shown in Figure IH i s the 
period -luminosity relationship for LMC Miras bv iFeast et alJ 
(1989) with the zero point as determined for galactic Miras by 
Whitelock & Feast (2000): 

Mk = -3.47 logP(days) + 0.85 (4) 
FigureEJshows several features: 
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Fig. 4. Distribution of LPVs in the Mk — logF plane, (a) 
Upper panel: stars are sorted by variable type. Filled circles: 
Miras. Open circles: SR/SRb. Crosses: SRa. (b) Lower panel: 
stars are sorted by chemical type. Filled symbols: carbon stars. 
Open symbols: oxygen stars. The dashed line is the P-L rela- 
tionship of Eq. 0. 

1 . To first order, the location of LPVs on the P-L diagrarn is a 
scatter plot, as shown also bvlBedding & Ziilstral(ll998h for 
oxygen Miras and bv lBergeat et alJ ( 200 ll) for'carbon stars. 

2. The Miras and SR are par tly-separated on this diagram (cf. 
also lBedding & Ziilstral l998). The Miras have longer peri- 
ods (no star with a period shorter than 120 days is a Mira), 
but both types of variables have similar luminosity ranges. 

3. The separation is not complete. Several SR are located in 
the Mira region, and vice versa. 

4. The pulsation characteristics of carbon and oxygen-rich 
Miras are indistinguishable. 

5. The SR P-L sequence, if there is one, has a different slope 
from that of the Miras. 

6. The partial overlapping of the Mira and SRb regions im- 
plies that some Miras may switch to SRb and vice-versa. 
As we discuss in the next section, there is ample evidence 
for mode-switching between Mira-Uke and SR-like pulsa- 
tion modes, and the incomplete segregation in Figure0]may 
be due to this. 

7. The (relatively small) number of SRa stars are spread ap- 
proximately equally between the Mira and SRb regions of 
the diagram. 



8. Both the Mira and SR variables have a maximum luminos- 
ity at K of Mk = — 8.2 ± 0.2. The maximum luminosity 
also holds for both oxygen and carbon stars, and is in agree- 
ment with the value found for LMC stars. 

3.1. The Period-Luminosity relationship 

We now use the data presented in Figurel^to calculate the mean 
P-L relation. We discuss the Mira and SRb samples separately 
and do not include the SRa. 

Many of the objects do not have significant parallaxes, even 
after reprocessing the Hipparcos data, and absolute magnitudes 
derived from the subset of objects in Tables ^ and |2l and de- 
picted in Figure |3 whic h are selected by par allax (and/or by 
vj/e(w)) will be biased (jLutz & KeUceJ 19731) . The correct es- 
timation of mean absolute magnitudes is then analogous to 
weighting the "detections" of parallax by the upper limits on 
the parallax on the assumption that the stars with detected and 
non-detected parallaxes are drawn from the same population. 
While the Hipparcos catalogue is magnitude limited and there- 
fore absolute magnitudes are also subject to Malmquist bias, 
the parallax for an LPV may become undetectable with respect 
to stars of similar brightness before the star does - a star with 
Mk = —8 at D = 1 kpc has V ^ 8 - and errors on the paral- 
laxes are therefore the major source of uncertainty. 

Although th e Astrometry Based Luminosity (ABL, 
lArenou & Luri 1999) is an unbiased estimator of the lu- 
minosity, it cannot be used directly to calculate absolute 
magnitudes which can then be used to derive a period- 
luminosity relation: the symmetry of the error bars on ABL 
does not propagate to the absolute magnitude. Instead of 
fitting the absolute magnitudes, we work in parallax space, 
where the errors are Gaussian. Therefore, instead of fitting 
Eq.©, we fit 



m = P-IO^IQ- 



(5) 



where m, the parallax in mas, can be used regardless of its 
sign. The parameters are fitted by chi-square minimization 
using non-linear optimization technique, and their uncer- 
tainties are obtained by mean of a standard Monte Carlo 
simulation. 

Figure |5l shows the resulting mean Mk — log P diagram 
for Miras. The data in that figure are strongly correlated: the 
linear correlation coefficient is 0.84, giving a ^ 99% chance 
of correlation. The least-squares fit to the P-L relationship for 
Miras gives 



Mk = -3.22(±0.38) logP(days) + 0.45(±2.40) 



(6) 



which is very consistent with that found by iFeast et all 
(1989) and is shown in Fig.|51 

Figure|5]also shows P-L data for three nearby OH/IR star s 
with reliable phase-lag distances dvan Langevelde et alJll990l) : 
WX Psc (IRCh-10011), RR Aql and R Aql. The last object is 
nearby enough to be measured by Hipparcos: the phase lag 
distance, 290 pc, is in reasonably good agreement with the 
Hipparcos parallax {zu^^ = 235 ± 52 pc). These data also 
show reasonable correspondence with the derived P-L slope. 
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log P (days) 

Fig. 5. Binned P-L data for Mira variables. The LMC relation- 
ship of tpeast et alj 11989) is shown by the dotted line, and the 
least-squares fit to the data by the solid line. Crosses: oxygen- 
rich Mira variables with distances from OH phase-lag measure- 
ments (see text). 
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The data for the SR variables also show a correlation, 
though less strongly: the correlation coefficient is 0.65, giving 
a >90% chance of correlation. The least-squares fit is: 

Mk = -1.34(±0.06) logF(days) - 4.5(±0.35) (7) 

The above res ults disagree witli witli tliose of 
[Bartlies et al.' ^1999') who divide their data into different 
samples, finding different relations for each sample, with 
slopes shallower than the LMC Mira slope. 

4. Variability and l\1ass Loss 

Miras with periods as long as nearly 2000 days have been 
found. These are all OH/IR stars with large mass loss rates. 
As discussed below, mass loss becomes very high for long- 
period variable stars, completely obscuring the star at visible 
and near-infrared wavelengths. Thus a large part of the period 
space occupied by LPVs is inaccessible to studies of the P-L re- 
lationship because parallax values for these stars are measured 
at present only at optical wavelengths. However, there are indi- 
cations that the P-L relationships so far found do not hold for 
very long period Miras. First, if it did, the predicted bolometric 
luminosity for a star with a 2000 d period would be 7 x 10"' L0, 
well in excess of the maximum AGB luminosity of 5 x lO'' Lf^, 
JPaczvhs ki 1970). Recently, 'He & Chen (2001) have examined 
the P-L relationship for OH/IR stars with phase-lag or kine- 
matic distances; there is a large amount of scatter in the P-L 
diagram, and the relationship may be different for objects in 
different IRAS spectral classes. Furthermore, the bolometric 



Fig. 6. Infrared color (S12/S2.2) versus period for oxygen stars 
(upper panel) and carbon stars (lower panel). Open symbols: 
Mira variables. Crosses: SR/SRb variables. The horizontal line 
shows the S12/S2.2 for a 2500 K black body. 

luminosities of many of these highly obscured, very long pe- 
riod stars are relatively modest, often < 10* L0. 

The CGCVS contains very few stars with periods longer 
than 400 day s. However , such stars ar e known, for example the 
OH/IR stars. 'Le Bertt^ (Il992l Il993l) measured near infrared 
light curves for several tens of oxygen-rich and carbon stars 
which were first discovered in infrared sky surveys. These ob- 
jects are heavily obscured, and almost all of them have periods 
longer than 400 days. Almost all of them vary with large am- 
plitudes and are therefore Mira variables. In Figure |6l the IR 
colors of these stars are compared with those from Tables^and 
2, showing the ratio of the 12/im/2.2/im flux densities versus 
the period (as noted in Section 2, this quantity is a good indi- 
cator of mass loss). The objects plotted are those from Tables 
[2and|2for which we have the appropriate m easurements, plus 
the stars for which iLe Bertr3 ( 1 19921 1 19931) measured periods 
(with the exception of the OH/IR stars, which are often at large 
distances and therefore have an unknown amount of interstellar 
extinction), an d, to round out the samples, the carbon star CIT 
6 [period from Alksnis '('1995')] and the oxygen supergiants VY 
CMa (Mira) and VX Sgr (SRb). 

Figure |6l shows that copious mass loss turns on more-or- 
less abruptly at periods of about 300 - 350 days for both oxygen 
and carbon stars. Further, the mass loss rate is much higher for 
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Fig. 7. Example COBE K band light curves for LPVs. T Cep 
and S Cep are Miras, R UMi is SRa, and the other stars 
are SR. A small number of data points with large variance 
have been removed. Both T cep and S Cep show the bump 
on the rising part of the h ght curve often seen for Miras 
jLockwood & Zwing|l9'7ll) 



Miras at a given period than it is for SR variables (the flattening 
in S12/S2.2 for oxygen/S stars at long periods is due to the sat- 
uration of the I2^m flux as the 9.7 fim silicate feature goes into 
absorption with increasing optical depth). Stars with high mass 
loss rates are essentially invisible at optical wavelengths. Thus 
the P-L relation is incomplete: stars with periods longer than 
about 400 days are absent. This part of period space will only 
become accessible by astrometry at the few microarcseconds 
level or better at infrared wavelengths. 



days: T Cep's period is 388 d, S Cep's 488 d. The partial COBE 
light curves are consistent with this period. 

Table |3] lists period data from DIRBE, where available, 
and from the literature, for stars with multiple periods, all 
of which are in the sam ple analyzed earlier in this paper. 
iKerschbaum et alJ (|2001) note that fewer than 50% of the pe- 
riods for SRb stars determined from their monitoring observa- 
tions agree with those from the CGCVS: the data in Table |3] 
support this conclusion. 

Table 3. Periods 



HIP 


GCVS 


var 


P 


Ref 


21479 


RDor 


SRb 


209 


DIRBE 








175,332 


Bedding & Ziilstra (19981 








338 


CGCVS 


28166 


BQ Ori 


SRa 


127,240 


Kiss et al. fl999) 








110 


CGCVS 


36288 


Y Lyn 


SRc 


133,205 


Szatmarv et al^ 








110,1400 


Percv et al. (2001 ) 








110 


CGCVS 


53085 


V Hya 


SRa 


530,6000 


KnaoD et al. C1999") 








529 


CGCVS 


63950 


FS Com 


SRb 


55, 600 


Percv et al. f200D 








58 


CGCVS 


65835 


R Hya 


M 


(314) 


DIRBE 








388 


CGCVS 


67419 


W Hya 


SRa 


274 


DIRBE 








382 


CGCVS 


68815 


6) Aps 


SRb 


121 


DIRBE 








119 


CGCVS 


70969 


YCen 


SRb 


(338) 


DIRBE 








180 


CGCVS 


80802 


RUMi 


SRa 


181 


DIRBE 








170,325 


Kiss et al. 0999) 








324 


CGCVS 


81747 


AXSco 


SR 


(39),(64) 


Kerschbaum et al. f2001) 








124,128 


Kerschbaum et al. f2001) 








138 


CGCVS 


81835 


S Dra 


SRb 


182 


DIRBE 








136 


CGCVS 


95173 


TSge 


SR 


(112) 


Kerschbaum et al. (2001) 








178,316 


Kerschbaum et al. f2001) 








166 


CGCVS 


100605 


UUDra 


SRb 


117 


DIRBE 








120 


CGCVS 


104451 TCep 


M 


112,>290 DIRBE 








388 


CGCVS 



5. Variability modes in LPVs 

Figure0illustrates some additional problems with studying the 
P-L relation for LPVs. It shows DIRBE 2.2/im light curves for 
six stars for which almost complete time coverage is available: 
two Miras (S Cep and T Cep), one SRa (R UMi) and three SRb 
(0 Aps, UU Dra and R Dor). The SR variables all have periods 
short enough that a complete phase information is available, 
and the SRa, R UMi, is indistinguishable in its variability char- 
acteristics from the other SR stars. Both Mira variables have 
periods that are longer than the COBE time coverage of 300 



Several analyses of very long time series data fo r individ- 
ual sta rs have recently appeared in the literature. Bedding et alJ 
(Il998l) find two periods, 175 d and 332 d, for R Dor, wifli 
switching between the two periods on a timescale of about 
1000 days. The DIRBE light curve has an intermediate pe- 
riod of 209 d, intermediate between the long and short periods. 
Bedding et al. ( 1998) note that the longer period, plus the star's 
well-determined absolute magnitude, place it on the LMC Mira 
P-L relation (note that it also lies on the Mira P-L relation found 
above) and suggest that the star switches pulsation modes be- 
tween SR and Mira variability, i.e. between third and first over- 
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tone radial pulsation, with Miras considered to be pulsating in 
the first overtone (Wood 1990; Willson 2000) - although see 
the discussion by Barthes ( 1998) . 

As pointed out by .Lebzelter & KissI (l200lh understand- 
ing LPVs requires both short- and long-term monitoring. 
Other studies of long tim e-series for LPVs include that of 
iHowarth & GreavesI ( 1200 ll) who find two periods for the Mira 
variable T Cas, 445 d and 222 d, whos e ratio suggest s that the 
longer period is the fundamental (cf. lBartheslfl998h : switch- 
ing between these two modes takes pla ce on a time scale of 
about 3000 d. Szatmar vet alJ(ll996l) and lBedding etal] ( ll998l) 
find that V Boo appears to have changed from a Mira to an 
SR variable. Icadmus et al.! also find evidence of mode 

switching for several stars. Sterken et al. ( 1999) show that the 
period of x Cyg is increasing. One star which does not ap- 
pear to undergo such behavior is V Hya, for which lKnapp et al] 
d 19991) find repeata ble and well establishe d periods of 530 d 
and 6000 d. Finally, iMarengo et al.l (I^OOl*) examine the struc- 
ture of circumstellar shells due to mass loss, in which the time- 
history of that mass loss may be seen, and suggest that essen- 
tially all LPVs undergo mode switching between the Mira and 
semi-regular phases. 

6. Discussion and Conciusions 

In this paper, we have examined the distribution of LPVs in 
nearby regions of the Galaxy in the P — Mk plane using dis- 
tances derived from Hipparcos parallaxes. We find P-L rela- 
tionships for both Mira and SRb variables. The slope for the 
SRb is much shallower than that for the Miras. The relation- 
ships found herein do not agree with some previous relation- 
ships found in the literature. The P-L relationship for Miras 
agrees well with that found in the LMC. 

We discuss information from the present paper and from 
the literature on period changes in LPVs and on mass loss and 
its dependence on the variability period. The evidence suggests 
that LPVs, at least those with periods < 400 d, change vari- 
ability modes back and forth between Mira and SR variation, 
with models suggesting that different radial pulsation modes 
become dominant. However, the P-L relationships found here 
do not extend to periods longer than about 700 days. 

If visible LPVs, i.e. those with low mass loss rates, switch 
between the Mira and SR modes, the shallowness of the P-L 
relationship for SRb relative to that for Miras is easy to under- 
stand: the period of a given star may decrease during the Mira 

SRb phase by a factor of 2-4 (depending on the pulsation 
modes) but the luminosity stays about the same. 

We also show that mass loss is strongly dependent on pe- 
riod and variability mode. Miras at a given period are losing 
much more mass than are SR, and the mass loss rate rises 
steadily with period for both Miras and SRb. Stars with peri- 
ods longer than about 400 - 500 days are losing so much mass 
that they become highly obscured at visible wavelengths and 
have no parallax information: thus a large part of period space 
is not available for studies of the P-L relations. These very long 
period variables are almost all Miras, suggesting that the mode 
switching characteristic of shorter-period red variables disap- 
pears as the period lengthens. 



Comparing the P-L relationships for LPVs from the litera- 
ture and this paper, it appears that the results are strongly de- 
pendent on the sample selection. The use of P-L relationships 
to determine distances to LPVs is therefore not likely to be 
useful in at least the near future. However, the maximum K 
absolute magnitude in the nearby (few 100 pc) region of the 
Galaxy is —8.2 ± 0.2, in good agreement with observations 
of the LMC. This corresponds to a bolometric luminosity of 
about 5000 - 6000 Lq, about the maximum possible luminos- 
ity for AGB stars. This result, in itself, shows the unlikeliness 
of a P-L relationship which extends to periods longer than 400 
- 500 days: stars with shorter periods are akeady close to the 
AGB luminosity limit. 

Further progress on studying the P-L characteristics of 
LPVs would be greatly aided by astrometric measurements at 
the few micro-arcsecond level in near-infrared bands, which 
would make possible the analysis of the entire period range. 
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